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United States Patent n9] [in 3,851,162 

Miinoz [45] Nov. 26, 1974 

157] ABSTRACT 

An input analog signal to be frequency analyzed is 
separated into N number of simultaneous analog sig- 
nal components each identical to the original but de- 
layed relative to the original by a successively larger 
time delay. The separated and delayed analog compo- 
nents are combined together in a suitable number of 
adders and attenuators in accordance with at least one 
component product of the continuous Fourier trans- 
form and analog signal matrices to separate the analog 
input signal into at least one of its continuous analog 
frequency components of bandwidth \/N times the 
bandwidth of the original input signal. Given the sepa- 
rated frequency components, the original analog input 
signal can be reconstituted by combining the separate 
analog frequency components in accordance with the 
component products of the continuous Fourier trans- 
form and analog frequency component matrices. The 
continuous Fourier transformation is useful for spec- 
trum analysis, filtering, transfer function synthesis, and 
communications. 


22 Claims, 10 Drawing Figures 
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CONTINUOUS FOURIER TRANSFORM METHOD 
AND APPARATUS 

The invention described herein was made by an em- 
ployee of the United States Government and may be 
manufactured and used by or for the Government for 
governmental purposes without the payment of any 
royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

The present invention relates in general to method 
and apparatus for performing Fourier transformations 
and more particularly to an improved continuous Fou- 
rier transformation and systems using same. 

DESCRIPTION OF THE PRIOR ART 

Heretofore, analog input signals have been analyzed 
by sampling the analog input signal at a multiplicity of 
time displaced intervals, converting such sampled data 
to discrete digital form in an analog to digital converter 
and performing a discrete Fourier transformation, typi- 
cally by the fast Fourier transform method, to convert 
the signal data in the time domain into its correspond- 
ing discrete frequency components in the frequency 
domain, thereby obtaining a discrete spectral represen- 
tation of the input signal to be analyzed. 

It is als o known from the prior art that on ce an a nalog 
signal has been separated into its corresponding dis- 
crete Fourier components in the frequency domain that 
a discrete representation of the original analog signal 
can be reconstituted by performing an inverse Fourier 
transformation on the discrete spectral Fourier compo- 
nents in a digital computer to derive a discrete repre- 
sentation of the original analog signal in the time do- 
main. The discrete fast Fourier transform method and 
apparatus is disclosed in an article “What is the Fast 
Fourier Transform?,” appearing in the I.E.E.E. Trans- 
actions on Audio and Electroacoustics, volume 
AU-15, No. 2, of June 1967, pages 45-55. 

One of the problems with the discrete Fourier trans- 
form is that it requires the use of a digital computer 
with associated analog-to-digital converters and the 
analysis performed by the computer comprises many 
complex multiplications, additions, and memory cy- 
cles. 

SUMMARY OF THE PRESENT INVENTION 

The principal object of the present invention is the 
provision of an improved Fourier transform method 
and apparatus and more particularly of a simplified 
continuous Fourier transform method which does not 
require a complex digital computer for implementa- 
tion.- 

In one feature of the present invention, an analog sig- 
nal to be analyzed is separated into N number of analog 
components each delayed relative to the original signal 
or the reference one of said signal components by a 
successively larger time delay, such reference and de- 
layed components are added together in accordance 
with at least one of the component products of the fast 
Fourier matrix (D matrix) and the delayed input signal 
component matrix F(i) to produce at least one contin- 
uous output function Cmil) corresponding to one of the 
frequency components of the analog input signal to be 
analyzed. 

In another feature of the present invention, the Fou- 
rier components in the frequency domain of an original 


2 

analog signal are combined in accordance with the 
component product of the discrete Fourier transform 
matrix (D matrix) and the Fourier components matrix 
(Cn(/) matrix) to reconstitute the original analog input 
5 signal. 

In another feature of the present invention, an analog 
input signal having a certain bandwidth is Fourier 
transformed into N number of Fourier frequency com- 
ponents each component having \IN X the bandwidth 
10 of the original input signal. Such separate Fourier fre- 
quency components are transmitted simultaneously 
over N number of different channels to a receiver. The 
received Fourier components are recombined in accor- 
dance with the inverse Fourier transform to reconsti-; 
15 tute the original input signal, whereby a relatively wide 
bandwidth signal is transmitted over N number of rela- 
tively narrow band individual channels. 

Other features and advantages of the present inven- 
tion will become apparent upon a perusal of the follow- 
20 ing specification taken in connection with the accom- 
panying drawings wherein: 

BRIEF DESCRIPTION OF THE DRAWINGS 

nG.‘ 1 is a signal flow graph depicting the continuous 
25 Fourier transformation of a complex (two-wire) analog 
input signal into four separate complex Fourier fre- 
quency components, 

RG. 2 is a signal flow graph depicting the continuous 
Fourier transformation of a real (one-wire) input signal 
30 to its respective Fourier frequency components, 

FIG. 3 is a physical implementation of the system of 
RG. 2, 

RG. 4 is an alternative physical implementation of 
the input portion of the circuit of RG. 3 and delineated 
35 by line 4 — 4, 

RG. 5 is a plot of the input signal to be Fourier ana- 
lyzed and the output Fourier components as a function 
of time and depicting the discrete Fourier transform 
output indicated by the dots superimposed upon the 
continuous Fourier transform output, 

FIG. 6 is a schematic circuit diagram for a physical 
implementation of a continuous Fourier transform cir- 
cuit similar to that of FIG. 3 except for eight time de- 
layed input comfjonents, 

FIG. 7 is a composite plot of band limited input signal 
data, frequency bands of the separate Fourier compo- 
nents, and the frequency response of the first Fourier 
frequency output coefficient Co, all as a function of fre- 
50 quency, 

RG. 8 is a schematic block diagram of a data trans- 
mission system incorporating features of the present 
invention, 

RG. 9 is a diagram similar to that of FIG. 8 for modi- 
fying an input signal, and 

RG. 10 is a schematic circuit diagram of an inverse 
continuous Fourier transform apparatus. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

The prior niethod of discrete Fourier transform will 
first be described followed by a description of an exten- 
sion thereof under the present invention to a continu- 
ous Fourier transform method and apparatus. In the 
discrete Fourier transform method, a band limited con- 
tinuous complex function of timefit) is considered to 
contain a maximum frequency component o)J2 
■ wherein w, is the rate at which the function /(t) is to be 
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sampled to produce a sampled data function J*{nT), 
where T = l/oj, seconds. A set of N sequentially or- 
dered magnitudes of/* (nT) where n = 0, 1, . . . , yv — 
1 , is selected and the discrete Fourier transform of this 
set is as follows: 

DFT (/«)) = c„.=i XJ 

n=0 

Eq. (1) 

where 

j= v=r 

m = 0, 1 , . . . A/-1 

and c„ are the N complex frequencies coefficients. 

Tlie inverse operation of computing the original data 
set from the coefficients is given as follows: 


fast Fourier transform is relatively efficient because 
using the fast Fourier transform 2N log N complex 
multiply operations plus the attendant addition and 
memory cycles are required, whereas a straight imple- 
5 mentation of the matrix equations (3) and (4) would 
require approximately complex multiply operations 
and the fast Fourier transform reduces this by a factor 
of 2 log N/N which is appreciable for large N. 

In the continuous transform method of the present 
10 invention, the discrete Fourier transform equations (3) 
and (4) are utilized and the input signal f(t) is consid- 
ered as a vector quantity having real and imaginary 
components. Because Equations (3) and (4) are inde- 
pendent of time, the independent variable F can be 
15 made a function of time F(t) to produce the following 
continuous Fourier transform: 

C(t)= 1/yV DFU). 

Eq. ( .S ) 


f* in 

m=0 


, 20 

Matrix Equation (5) may be expanded and written as 
Eq. (2) follows: 


Co(0 

C,{() 

C,(») 


= \IN 


I I 

I e> 

1 fV 


1 

e' 


I 

gis-iti 


un 

f,U) 

Ml) 


Eq. (6) 
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where n =0, 1, . . . , N—l 
These equations can be written in matrix form in the 
following way: 


and 


C= \/N D F 


F=DC 


Eq. (3) 


Eq. (4) 


30 where 

As in the fast Fourier transform a considerable sim- 
plification in the computation is achieved if N is se- 
35 lected to be a binary number, i.e., a power of 2 such as 
2,4,8,16, etc. This simplifies F and greatly reduces the 
complexity. The matrix Equation (6), when solved for 
the complex frequency coefficients €„([), produces a 
set of equations as follows: 


C„(0=l//Vl/o(0+/,(() +/,(/) 

C,U)=UN\Mi)+e>Ml) -Ml) 

C,(()=IW|/,(()-/,(r) +MO 


■•• /«-,(()! (7) 

( 8 ) 

(01 (9) 

(10) 


where C — 


and F — 


fo 

C'l 


1 — 

1 

1 

. i 


D = 

. . ^ 

1 gifi- 1 b 

giis-iti . 



f> (o) 

/* (T) 


f ((/v-n 7) 


50 


For N = 4, these equations reduced to: 
Co(r) = \/N [foiO +/,(/) +Mt) +Mt)] 

Qio = \/N Ifoit) + e>Mt) -MO - e>Mt)] 

Q(r) = MN \foit) -MO +MD -MD] 


Eq. ( II ) 
Eq. (12) 
Eq. (13) 


The discrete fast Fourier transform is obtained for 
Equations (3) and (4) by recognizing and using the 
properties of the D matrix to simplify the computa- 
tions. Since all entries in the D matrix are periodic 
roots of unity, and symmetric about the main diagonal, 
a great simplification is possible where is a power of 
2. The highest degree of symmetry and hence the great- 
est simplification is permitted when N fulfills this re- 
quirement. The matrix equations of (3) and (4) are 
typically solved with a digital computer programmed to 
perform the fast Fourier transform. The program in- 
volves many complex multiplications and additions. Al- 
though the computations are relatively complex the 


C,(r) = MN [/o(r) - e'/,(r) -/(O + e^/Cr)] Eq. (i4) 

The signal flow graph for Equations 1 1-14 is illus- 
trated in FIG. 1 . In the signal flow graph, eaeh node 
represents a variable, and the arrows terminating at 
^0 that node originate at the nodes whose variables con- 
tribute to the value of the variable at that node. The 
contributions are additive, and the weight of each con- 
tribution, if other than unity, is indicated by the con- 
stant written close to the arrowhead of the transmis- 
^5 sion. Equations 1 l-i4 are separated into real and imag- 
inary components as follows: 
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ReQ,(t) = \!N {ReW) + RefM + Refi{t) + Ref^im 
- -- ■ (15) 

/mC„(0 = 1//V (/m/o(f) + ImfM + ImfM + /mMO] 

(16) 

ReQU) = l/N IReMO - lmf,(0 - Re/.O) + /m/,(0] 5 

(17) 

/mC,0) = 1/N l/m/oO) + Ref,(t) - ImUt) - ReMi)] 

(18) 

ReQU) = 1/A/ [RefoU) - RefM + RefM - /?£■/,(/)] 

( 19 ) 10 

ImQil) = UN [ImfaU) ~ Imft(t) + ImfiU) — lnif,(t)] 

(20) 

ReQO) = UN [RefoU) + - RefM ~ ImhU)] 

I (21) 

; ImQU) = \!N [/m/o(/) - RefM - hnUU) + RefzU)] 15 

, ( 22 ) 

Equations 15-22 could be implemented with hard- 
ware, however, considering a one dimensional input; 
signal having only a real and no imaginary component,: 
Equations 1 5-22 are greatly simplified by recognizing, 20 
that: 

tmfoU) = ImfiU) = ImfiU) = Imfsit) = 0 ‘ 

Thus: 

ReCoO) = 1/A/ [ReUt) + Re/M + RefM + RefAl)\ ^5 
/mC„(M = 0 

Kq. (241 

ReCM - ReC^U) = 1/A/ [Ref^U) ~ ReMO] Eq. (25) ^0 
ImQU) = -ImCzU) = UN [Refi(t) - /?e/3(0] Eq. (26) 

ReCiU) = 1//V [/?c/o(f) - RefM -1- Ref^U) ~ ReMO) 35 

Eq. (27) 


ImCtU) = 0 


Eq. (28) 


40 


ReC 3 {t) - UN [RefoU) ~ Ref 2 U)] Eq. (29) 
ImQU) = UN [- RefU) + ReMO] 

Eq. (30) 


45 


50 


55 


Since ImCoU) and ImC 2 (t) = zero, no hardware imple- 
mentation for Fourier frequency components is neces- 
sary or possible. Therefore, equations 23-30 reduced 
to only four independent equations, namely Equations 
1 3, 25, 26 and 27, the signal flow graph for which is de- 
picted in FIG. 2. In the continuous transform, /o(t) is 
identical to the reference input signal /(t) and /,(t) is 
identified with foU+T),Mi) is identified with foU-UlT) 
and foU) is identified with /o(/-t-37^. 

The signal flow graph of FIG. 2 is implemented with 
hardware as shown in FIG. 3 to provide a continuous 
Fourier transform apparatus. More particularly, the 
input function of time/(r) (magnitude) is first divided 
by N in attenuator 1 1 to provide a first reference input 
signal foU)!N. UN is a scaling constant equal to the re- 
ciprocal of the number of Fourier frequency compo- 
nents. A sample of the output of the attenuator 1 1 is 
fed to a first delay 12 to be delayed by a time T. The 
output of the delay 12 forms the second input signal, 
namely fU)!N. A sample of the output of (ielay 12 is . 
fed to a second delay 13 which delays the signal passing 
therethrough by the same time T such that the output 


of the second delay 13 forms the third input signal, 
namely ^(0/ A/. An output of the second delay 13 is fed 
to a third delay 14 for delaying the fourth signal by a 
time T relative to the third signal to produce the fourth 
input reference signals namely fsiON. The continuous 
input signals are combined in the inverting operational 
amplifiers 15 in accordance with the signal flow graph 
of FIG. 2 to produce the four outputs, namely, 
ReCoU). ReCiU). ReC,(t) = ReC,U). ImCfr) = 

— lmC:i[t) . - — 

One relatively simple apparatus for deriving the time 
delayed input signals for the input to the circuit of FIG. 
3 is as shown in FIG. 4. More particularly, the input sig- 
nal _/(/) to be Fourier analyzed is recorded on a mag- 
netic recording tape 16. The magnetic recording tape 
16 is drawn from a supply reel 17 onto a take-up reel 
18 across a number of pick-up heads 19 each head 19 
being equally spaced from the adjacent head in the di- 
rection of movement of the recording tape 16. The out- 
puts from each of the downstream hea(Is 19 is succes- 
sively delayed by a successively larger time compared 
to the output of the upstream head or reference head. 
The respective output signals are indicated above each 
of the respective heads. 

Referring now to FIG. 5, there is shown the signal to 
be analyzed yiO and its Fourier frequency components 
for A/ = 4 and corresponding to the four outputs of the 
circuit of FIG. 3. In FIG. S,fU) is shown being sepa- 
rated into A/ = 4 components and sampled four times 
per time window where the time window No. 1 is de- 
fined as 0 to 4T. The information contained in the four 
components ReCo(t), ReQit), ReCft) and ImCiU) is 
no greater nor less than that contained in the four origi- 
nal samples because the continuous Fourier transform 
process is conservative, that is, it permits information 
neither to be created nor destroyed. Since only one 
unique set of Fourier frequency components per win- 
dow is allowed to characterize fit), according to the 
sampling theorem, each Fourier frequency component 
must then contain only 1 IN X the bandwidth limit of fit 
) or a maximum frequency content of o,l2N = u)J8 Hz. 
This accounts for the bandwidth reduction in the dis- 
crete or sampled data case and this has a direct and im- 
mediate implication for the continuous transform 
mechanism. 

The continuous transform method is related to the 
series of discrete Fourier transforms in successive time 
windows 1, 2, ... p and at the limit of this process, i.e., 
when p goes to infinity and T goes to zero, the band- 
width relationships are invariant for both p and T. 
Therefore, they obtain in the limit and apply equally 
well for the continuous transform case, i.e., the band- 
width of any Fourier component Qit)= UN bandlimit 
of /(f). 

TTie continuous Fourier transform mechanism of the 
present invention is suited to many practical applica- 
tions including filtering, spectral analysis, and transfer 
function modification in the frequency or time domain. 
In short, it can do all that the theory permits of any 
Fourier transform mechanism and for some applica- 
tions it results in very simple and inexpensive hardware. 

Referring now to FIG. 6 there is shown a Fourier 
transform circuit similar to that of FIG. 3 where N = 8. 
As can be seen in FIG. 6, the hardware implementation 
involves merely attenuators and amplifiers as con- 
tra-sted with the prior art discrete Fourier transform 
mechanism that requires complex multiplication, addi- 



3,851,162 


7 

tion and memory cycles. The frequency characteristics 
of the continuous transform mechanism for the circuit 
of FIG. 6 is shown in FIG. 7 where N = 8 and lmf(t) 

= 0. The information band ofJ{t) from zero to oijl, is 
separated into five equibandwidth parts by the output 5 
frequency components Co through as indicated by 
frequency bands 0 through 4. This operation of produc- 
ing Fourier frequency components on a continuous 
basis is identical to the operation of filtering with the 
added advantage that an exact measurement of the 10 
phase of the incoming signal with respect to oi, is possi- 
ble for each band. It can be seen that the Co component 
acts as a low pass filter and Ci acts as a bandpass filter, 
etc., each band having a bandwidth of (+)o)J2N 
around their respective center frequencies of mojJN 15 
where m = 0, 1 . . . . (/V— 1). All this filtering action 
takes place simultaneously and continuously and there- 
fore constitutes a parallel set of A//2 filters or a spec- 
trum analyzer. If only one filter is needed, then the 
hardware can be simplified greatly. 20 

At frequencies above w,/2 the action of filtering, un-‘ 
like that of ordinary filters, repeats itself in a cyclical 
manner as indicated by the frequency band containing 
the same number designation in FIG. 7. This was the 
reason for band limiting the input signal f{t) to avoid 25 
signal aliasing and the consequent confusion of data 


8 

similar to that of FIG. 8 with the exception that variable 
attenuators 24 are placed in each of the respective 
channels between the continuous Fourier transform 
mechanism 21 and the inverse continuous Fourier 
transform mechanism 22 for modifying the frequency 
components in the frequency domain. The reconsti- 
tuted modified output signal G(t) appears at the output 
across a load resistor 25. Convolutions in frequency or 
time can be performed by using two continuous trans- 
form systems of the type shown in FIG. 9 and multiply- 
ing the input signals or the frequency components re- 
spectively and reconstituting or transforming these 
products using another transform mechanism. 

Once an input signal in the time domain has been 
separated into its continuous fast Fourier transform 
components Cm(/) in the frequency domain these com- 
ponents can be recombined to form the original signal 
by means of the inverse Fourier transform, the dis- 
crete version of which is shown in Eq. (4). In the con- 
tinuous Fourier transform. Equation (4) is rewritten as; 

F(t) = DCU) 

Eq. (.11 ) 

Equation (31) is solved for f,U) = AO and for A = 4 
as follows; 



Eq. (32) 


that results. However, for some applications this action 
may be desirable for use in a periodic filter as, for ex- 
ample, in the magnetic amplifier or flux gate magne- 
tometer where the output signals are contained in all 
even harmonics of the drive signal and it is necessary 
to reject the odd harmonics and extraneous noises in 
other bands. 

Many types of periodic filters are possible such as 
even and odd harmonics as well as third, fifth, eleventh, 
thirteenth, nineteenth . . . harmonic types, each se- 
lected as combinations of frequency bands containing 
in the same band number or numbers. This action can, 
of course, be accomplished simultaneously for different 
combinations and may be useful for human speech and 
music signal analysis. 

Returning now to band limited information input sig- 
nals, the continuous transform mechanism can be used 
to make a 'communication system capable of sending 
high bandwidth information over a plurality of N low 
bandwidth lines or channels. Such a transmission sys- 
tem is shown in FIG. 8. In FIG. 8, the input signal y(/) 
is continuously Fourier transformed in Fourier trans- 
former 21 in a manner as previously described with re- 
gard to FIG. 3. The output of transform mechanism 21 
is N frequency components with the attendent band- 
width reduction of 1/A for each of such Fourier fre- 
quency components. The components are then trans- 
mitted over separate channels 1 through A to be recon- 
stituted at the receiving end via an inverse continuous 
transform mechanism 22 into a facsimile of the original 
signal. 

Referring now to FIG. 9 there is shown an apparatus 
23 adapted for transfer function modification in the 
frequency domain. More particularly, the system is 


From Equation (32); 

AO = Co(0 + C,(0 + Q(4) -f Q(0 Eq. (33) 

Taking into account real and imaginary components, 

AO — ReCo(t) + ImCoiO + ReCi(t) + ImCi(t) -f 
ReC 2 (t) + ImRtU) + ReQiO + ImRaiO Eq. (34) 

Assuming the imaginary part of F{0 =0 then, ImC^U) 
= ImCziO = 0 and ReCiiO = ReCsiO and ImQ(0 = 
—ImC^iO then; 

AO = ReCoiO + 2 ReC,(/) + ReC^iO Eq. (3.s) 

Equation (35) is implemented by the hardware of FIG. 
10 and consists only of the positive sum of three com- 
ponents. The amplifier 22 is a non-inverting opera- 
tional amplifier with four input terminals each of a -f 
unity gain. The input signal ReCi(0 is connected into 
two input terminals to yield a gain of two. The output 
of the inverse continuous transformer 22 is the recon- 
stituted original input signal AO- 

The continuous Fourier transform mechanism of the 
present invention is useful for filtering, signal process- 
ing and communications. In filtering, it may be utilized 
for lowpass filtering, bandpass filtering, parallel filter- 
ing and chorus filtering (a selected set of harmonically 
related frequencies and combinations), specialized fil- 
tering of spectral signatures, and phase lock filtering or 
phase discrimination such as that used in space com- 
munication, and it is also useful for complex signal fil- 
tering. 

The advantage of the continuous Fourier transform 
apparatus of the present invention as applied to filter- 
ing is that very accurate frequency selection is obtained 
there being no practical limit on the minimum band- 
width. The filtering circuits are inexp>ensive to imple- 


35 


40 


45 


50 
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ment in hardware as the hardware comprises merely 
amplifiers and in some cases attenuators. 

In signal processing, such as in a system described 
with regard to FIG. 9, the continuous Fourier transform 
apparatus of the present invention is useful for transfer 5 
function modification in the frequency domain. It may 
also be utilized for spectrum analysis of continuous 
data (quasi stationary and non stationary time series 
produce deterministic results). It may also be utilized 
for serial-to-parallel data conversion and parallel-to- 10 
serial data conversion. 

In communications, such as the system described 
with regard to FIG. 8, the continuous transform mecha- 
nism is useful for bandwidth reduction in the individual 
communication channels. It may also be useful for 15 
video tape recording utilizing multiple tape heads. 

What is claimed is; 

1. In a method for continuous analog analysis for a 
time varying input function /(t) in the time domain into 
at least one of its Fourier Frequency components Cm(i) 
the steps of: 

continuously producing A/ signal components fJit)IN 
ranging from/o(/)/A^ X.ofs-x{t)!N wherein each sig- 
nal component except /o(/)/Af is delayed an amount 
T from the previous signal component; 

continuously combining separated ones of said signal 
components L/lnCO] fmWIN together so as to pro- 
duce at least one continuous output function C„(r) 
corresponding to one of the Fourier frequency 
components of the time varying input function^/) 30 
to be frequency analyzed. 

2. The method of claim 1 wherein the step of combin- 
ing said separated components fmU)IN together so as 
to produce at least one output function CmU) corre- 
sponding to one of the Fourier frequency components 35 
of the time varying input signal /(/) comprises the step 
of, combining the time delayed continuous signal com- 
ponents /m(t)//V together according to at least one of 
the component products of the matrix product: DF(t) 
where D is the discrete Fourier transform matrix and 40 
F(0 is the delayed input signal component matrix. 

3. The method of claim 2 wherein the step of combin- 
ing the delayed signal components /m(t)//V according to 
at least one component product of the matrix product 
includes the step of, adding the delayed signal compo- 45 
nents fmiOIN according to at least one of the compo- 
nent products of said matrix product. 

4. The method of claim 2 wherein D is defined by the 
matrix: 


jit) = DC(t) 

where D is the discrete Fourier transform matrix and 
C(r) is the matrix for the continuous separate Fourier 
frequency components. 

7. The method of claim 6 wherein the step of combin- 
ing a plurality of the continuous time varying Fourier 
frequency components Cm(r) together so as to reconsti- 
tute the original time varying function fit) includes, 
combining the separate Fourier frequency components 
C„it) according to the product of the matrix Equation: 

Jit) = DC{t) 

where D is the discrete Fourier transform matrix and 
Cit) is the matrix for the continuous separate Fourier 
frequency components. 

8. The method of claim 7 wherein the D matrix is de- 
fined by the matrix equation; 
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9. The method of claim 8 where 


C(0 = 


Cod) 
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5. The method of claim 1 wherein W is a power of 2. 

6. The method for reconstituting a time varying func- 
tion /(r) in the time domain from a plurality of N num- 
ber of its continuously time varying CmU) Fourier fre- 
quency components in the frequency domain the steps 
of: 

continuously combining the separate Fourier fre- 
quency components Cmit) according to the prod- 
uct of the matrix equation 
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10. In an analog apparatus for continuous frequency 
analysis of a time varying input function /(f) in the time 
domain into at least one of its Fourier frequency com- 
ponents Cm(t): 

means for continuously generating N signal compo- 
nents 

fmiOlf^ ranging from foU)/N to /,_,(»)/ A/ wherein each 
signal component except /o(r)/A/ is delayed an amount 
T from the previous signal component; 

means for continuously combining said signal com- 
ponents fmiOIN together so as to produce at least 
one continuous output function C„it) correspond- 
ing to one of the Fourier frequency components of 
the time varying input function /(f) to be frequency 
analyzed. 

11. The analog apparatus of claim 10 wherein the 
means for combining said separated signal compionents 
fmU)IN together so as to produce at least one output 
function CmU) correspionding to one of the Fourier fre- 
quency components of the time varying input signal 
jit) includes, means for combining the time delayed 
continuous signal components fm(t)IN according to at 
least one of the component products of the matrix 
product: 

DFU) 

where D is the discrete Fourier transform matrix and 
Fit) is the separated input signal component matrix. 

12. The analog apparatus of claim 11 wherein said 
means for combining the delayed signal components 
according to at least one component product of the ma- 
trix product includes, means for adding the delayed sig- 
nal components according to at least one of the compo- 
nent products of said matrix product. 
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13. The analog apparatus of claim 11 where D is de- 
fined by the matrix: 
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14. The analog apparatus of claim 10 wherein A/ is a 
power of 2. 

15. Apparatus comprising: 

first analog means for continuously producing A/ time 
varying Fourier frequency components Cm(r) in the 
frequency domain of a time varying function /(.r) in 
the time domain; 

said first analog means comprising: means for contin- 
uously generating N signal components fm(t)IN 
ranging from fa{t)!N to fs-i(t)!N wherein each sig- 
nal component except fo(t)IN is delayed an amount 
T from the previous signal component; and means 
for continuously combining said signal components 
fm(t)IN together so as to produce N varying Fou- 
rier frequency components C’m(t); 

channel means for conveying each time varying Fou- 
rier frequency component to an inverse Fourier 
transformer, said transformer comprising second 
analog means for combining the continuously time 
varying Fourier components C„(t) together so as to 
reconstitute the original time varying function 

16. The apparatus of claim 15 wherein said means for 
combining a plurality of the continuously time varying 
frequency components C„(/) together so as to reconsti- 
tute the original time varying function fit) includes, 
means for combining the separate frequency compo- 
nents Cm(t) according to the component products of 
the matrix product 

DCO) 

where D is the discrete Fourier transform matrix and 
C{t) is the matrix for the continuous separate Fourier 
frequency components. 

17. The apparatus of claim 16 wherein: 
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19. In an analog apparatus for continuous frequency 
analysis of a time varying input function^/) in the time 
domain into at least one of its Fourier frequency com- 
ponents Cm{i): 

5 means for continuously generating N signal compo- 
nents fmOVN ranging from foU)IN to /v_,(/)/A/ 
wherein each signal component except fo(t)/N is 
delayed an amount T from the previous signal com- 
ponent; 

10 said generating means including a tape recorder with 
a recording head and (N—2) staggered playback 
heads, said signal component fo{t)/N being con- 
nected to said recording head and each of the re- 
maining signal components being produced at one 
15 of said playback heads; 

means for continuously combining said signal com- 
ponents fmiOIN together so as to produce at least 
one continuous output function Cm(0 correspond- 
ing to one of the Fourier frequency components of 
20 the time varying input function^/) to be frequency 
analyzed. 

20. An analog apparatus as set forth in claim 19 
wherein said means for continuously combining signal 
components consists of a network of operational ampli- 

25 fiers. 

21. Apparatus comprising: 

first analog means for continuously producing N time 
varying Fourier frequency components Cm( t ) in the 
frequency domain of a time varying function in 
the time domain; 

said first analog means comprising: means for contin- 
uously generating N signal components fmU)IN 
ranging from fo(t)/N to fs-i(t)/N wherein each sig- 
nal except fo(t)IN is delayed an amount 7' from the 
previous signal component; and means for continu- 
ously combining said signal components fm(t)IN 
together so as to produce N varying Fourier fre- 
quency components Cm(/); 
said generating means including a tape recorder with 
a recording head and (N—2) staggered playback 
heads, said signal component fo(t)/N being con- 
nected to said recording head and each of the re- 
maining signal components being produced at one 
of said playback heads; 

channel means for conveying each time varying Fou- 
rier frequency component to an inverse Fourier 
transformer, said transformer comprising second 
analog means for combining the continuously time 
varying Fourier components Cm(t) together so as to 
reconstitute the original time varying function /(t); 


30 


35 


40 


45 


50 


18. The apparatus of claim 17 wherein: 


C(/) = 


c „(0 

c,(t) 


Cv-i(0 


said combining means consisting of a network of op- 
erational amplifiers. 

22. Apparatus as set forth in claim 21 wherein said 
55 inverse Fourier transformer consists of a network of 
operational amplifiers and attenuators. 
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